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Abstract 

I review recent progress in defining probability distributions in the inflationary multiverse. 
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I. INTRODUCTION 



It has long been a dream of particle physicists to derive the values of all constants of 
Nature from a fundamental theory. With the development of string theory in the last few 
decades, it seemed for a while that we were getting closer to that goal. String theory is 
our best candidate for the fundamental theory, and there has been great enthusiasm and 
hope that it will yield a unique set of constants. This hope, however, appears to have been 
dashed. It now appears that string theory has a multitude of solutions describing vacua with 
different values of the low-energy constants. The number of vacua in this vast "landscape" 
of possibilities can be as large as 10 

In the cosmo logical context, high-energy vacua will drive exponential inflationary ex- 
pansion of the universe. Transitions between different vacua will occur through quantum 
tunneling, with bubbles of different vacua nucleating and expanding in the never-ending 
process of eternal inflation. As a result, the entire landscape of vacua will be explored. 

If indeed this kind of picture describes our universe, then we will never be able to cal- 
culate all constants of Nature from first principles. At best we may only be able to make 
statistical predictions. The key problem is then to calculate the probabihty distribution for 
the constants. It is often referred to as the measure problem. 

n 

The probability Pj of observing vacuum j can be expressed as a product |^ 

P^. = P^Wnf^), (1) 

where Pj^^ is the fraction of volume occupied by vacuum of type j and nj"^*"* is the number 
of observers per unit volume. The distribution then gives the probability for a randomly 
picked observer to be in a vacuum of type j. 

The density of observers n^"'^'^^ cannot at present be calculated, but in many interesting 
cases it seems reasonable to approximate it as being proportional to the density of suitable 
stars, which is in turn proportional to the fraction of matter clustered in giant galaxies 

BflQ. 

The volume factor Pj^^ presents a problem of a very different kind: the result depends 
very sensitively on the choice of a spacelike hypersurface (a constant-time surface) on which 
the distribution is to be evaluated. This problem was uncovered by Andrei Linde and his 

collaborators when they first attempted to calculate volume distributions 0,y, Q. It 
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FIG. 1: A schematic conformal diagram for a comoving region in an eternally inflating universe. 
Bubbles of different vacua are represented by different shades of gray. The upper boundary of the 
diagram is the future timelike infinity. A surface of constant global time S cuts through the 
entire region and intersects many bubbles. 

eluded resolution for more than a decade, but recently there have been some promising 
developments, and I believe we are getting close to completely solving the problem. The 
purpose of this paper is to review the new proposals for Pj^^ ■ 

II. PROBLEM WITH GLOBAL-TIME MEASURE 

The spacetime structure of an eternally inflating universe is schematically illustrated in 
Fig.l. The bubbles expand rapidly approaching the speed of light, so their worldsheets are 
well approximated by light cones. If the vacuum inside a bubble has positive energy density, 
it becomes a site of further bubble nucleation; we call such vacua "recyclable". Negative- 
energy vacua, on the other hand, quickly develop curvature singularities; we shall call them 
"terminal vacua". 

The diagram represents a comoving region, which is initially comparable to the horizon. 
The initial moment is a spacelike hypersurface Sq, represented by the lower horizontal 
boundary of the diagram, while the upper boundary represents future infinity, when the 
region and all the bubbles become infinitely large. How can we find the fraction of volume 
occupied by different vacua? A natural thing to do is to consider a spacelike hypersurface S, 
which cuts through the entire region, as shown in the figure. If t is a globally defined time 
coordinate, then all surfaces t = const will have this property. One can use, for example, the 
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proper time along the "comoving" geodesies orthogonal to the surface Sq.^ Alternatively, one 
could use the so-called scale factor time, defined as a logarithm of the expansion factor along 
the comoving geodesies, or any other suitable time coordinate. Once the time coordinate 
is specified, one can find the fraction of volume occupied by different vacua on the surface 
t = const and then take the limit t ^ oo. 

Unfortunately, as I have already mentioned, the result of this calculation is sensitively 
dependent o„ one. choice of the t.. coordinate g The .ea.on . that the volume ot an 
eternally inflating universe is growing exponentially with time. The volumes of regions filled 
with all possible vacua are growing exponentially as well. At any time, a substantial part 
of the total volume is in new bubbles which have just nucleated. Which of these bubbles 
are cut by the surface depends on how the surface is drawn; hence the gauge- dependence 
of the result. Since time is an arbitrary label in General Relativity, none of the possible 
choices of the global time coordinate a ppe ars to be preferred. For more discussion of this 



gauge-dependence problem, see 
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III. A POCKET-BASED MEASURE 



It is now becoming increasingly clear that the solution to the problem lies in the direction 
of using the local definition of time within individual bubbles, or, as Alan Guth called them, 
"pocket universes" . I will first discuss how this works in the simplest case, when there is a 
single type of bubble; the general situation will be considered in the next section. 

Suppose we have an eternally inflating universe filled with a metastable false vacuum 
F, which decays to the true vacuum T through bubble nucleation. The vacuum energy 
is thermalized within the bubbles, and in due course observers evolve there and measure 
the constants of Nature. Suppose further that there is a scalar field X, which affects the 
values of some constants and has a very slowly varying potential U{X). The values of X 
are randomized by quantum fluctuations during inflation, so X is slowly varying in space 
within each bubble. We shall assume that the slope of U (X) is so small that time variation 
of X is negligible during the epoch when the observers are present. Our goal is to find the 

^ The term "comoving" is used very loosely here, since the vacuum does not define any rest frame. Any 
congruence of geodesies orthogonal to a more or less flat spacelike surface Sq can be regarded as "comov- 
ing". 
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FIG. 2: A comoving region with a single type of bubble. Bubbles differ by details of the scalar 
field distribution, but are statistically identical. Each pocket is internally an infinite open universe. 
Constant-time surfaces within pockets, shown by dashed lines, are infinite spacelike surfaces. 



distribution 

P^^\X)dX, (2) 

which gives the volume fraction occupied by regions where the field is in the interval between 
X and X + dX. 

A measure based on a global time coordinate runs into the same problem as in the case of 
a discrete distribution pj^-*: the result for P^^UX) is gauge- dependent. There is, however, 

nnn n 

a simple way around this difficulty [1^ • It is well known [12l| that bubble interiors 

appear to their inhabitants as self-contained infinite universes of negative curvature. A 
natural definition of the time coordinate in such a universe is to identify it with one of the 
physical variables, e.g., the energy density. Each pocket universe will then have its own set 
of infinite constant-time surfaces (see Fig. 2). 

The p..„po.a. of QHQ is to calculate the d.t.but.o. w.Mn a s.g.e poCet 
universe. It does not matter which one, since all pocket universes are statistically equivalent: 
they all have the same volume distribution of X. 

Once a pocket universe is selected, we choose some constant-time surface within it as our 
reference surface S. The next step is to find the volume distribution of X on S. Some care is 
needed here, since S is an infinite hypersurface. If one calculates the distribution in a region 
of finite size and then takes the limit of size going to infinity, the result may depend on the 
limiting procedure. To avoid this danger, one has to make sure that the shape of the region 
is not correlated with the distribution of X. The simplest choice is to use a "spherical" 
region, defined as a set of points with a geodesic distance r < R from a given center, with 



the limit R oo taken afterwards 



Since we assume that X does not change in time, its value remains constant along the 
comoving geodesies, and thus the distribution Q is time-independent if we think of it 
as a comoving- volume distribution.^ The density of observers n^^^^^X) should then also 
be understood as the number of observers per unit comoving volume. If we include all 
present, past and future observers, then n^°'''^\X) is also time- independent. Moreover, the 
full distribution 

dP{X) oc P(^) (X)n(°^^) (X) (3) 

does not depend on the choice of the reference hypersurface S, as long as P^^\X) and 
n^°^^\X) are evaluated on the same hypersurface. 

Some analytic and numerical techniques for calculating P^^\X) have been suggested in 
3 H Q. The=e techniques a.e not yet Mly developed, and there a.e some mte.estmg 
issues that still need to be addressed (e.g., the "ergodic conjecture" in Ref. PJ]). But as a 
matter of principle, the problem appears to have been solved in the case of identical pockets. 



IV. COUNTING POCKETS 



If there are many different types of pockets, then it is clearly not sufficient to consider 
a single bubble. We have to learn how to compare the numbers of observers in different 
bubbles. Since the bubbles are disconnected from one another, we have to define a comoving 
length scale Rj on which observers are to be counted in bubbles of type j. In addition, some 
bubbles may be more abundant than others, and we have to introduce a frequency factor 
Pj characterizing the relative abundance of different bubbles. The full expression for the 

^ An alternative definition of the distance is in terms of the area element, dA = r^dfl, with solid angle dfl 
defining a bunch of radial geodesies. It is possible that different definitions of r may yield different volume 
distributions for X, even in the limit r oo. The reason is that in a space of negative curvature, most of 
the volume of a sphere is near its surface (f is an exponentially growing function of r). The distributions 
for X may differ if the spacetime geometry is influenced by the local value of X. (This would be analogous 
to the gauge-dependence of a global-time measure.) This problem does not arise if the comoving sphere is 
defined at an early time, when the open FRW universe inside the bubble is dominated by curvature (see 
Sec. IV.C). 

^ The physical- volume distribution will generally evolve, since the local value of X may affect the expansion 
rate of the universe. 
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volume distribution is then given by [19| 

« p^R]- (4) 

Here I assume for simplicity that there are no continuous fields X that can vary within 
bubbles. The discussion can be easily extended to include such variables (see Sec. IV. D). 
The question now is: How do we define pj and Rjl We shall consider them in turn. 



A. Bubble abundance pj 

The definition of pj is a tricky business, because the total number of bubbles is infinite, 
even in a region of a finite comoving size.^ We thus need to introduce some sort of a cutoff. 
Here I shall review the procedure recently proposed in which has some attractive 
properties. 

The proposal is very simple: count only bubbles greater than a certain comoving size e, 
and then take the limit e ^ 0. To define the comoving size, one has to specify a congruence 
of "comoving" geodesies emanating from some initial spacelike hypersurface Sq. As they 
extend to the future, the geodesies will generally cross a number of bubbles before ending 
up in one of the terminal bubbles, where infiation comes to an end. There will also be 
a (measure zero) set of geodesies which never hit terminal bubbles. The startin g p oints 



21 



of these geodesies on So provide a mapping of the eternally infiating fractal 
consisting of points on where infiation never ends. In the same manner, each bubble 
encountered by the geodesies will also be mapped on Eq, and we can define the "size" of a 
bubble as the volume of its image on Sq. (The volume of a bubble is calculated including 
all the daughter bubbles that nucleate within it.) 

In an infiating spacetime, geodesies are rapidly diverging, so bubbles formed at later 
times have a smaller comoving size. (The comoving size of a bubble is set by the horizon 
at the time of bubble nucleation.) The bubble counting can be done in an arbitrarily small 



^ The problem of calculating pj is somewhat similar to the question of what fraction of all natural numbers 
are odd. The answer depends on how the numbers are ordered. With the standard ordering, 1, 2, 3, 4, 
the fraction of odd numbers in a long stretch of the sequence is 1/2, but if one uses an alternative ordering 
1, 3, 2, 5, 7, 4, the result would be 2/3. One could argue that, in the case of integers, the standard 
ordering is more natural, so the correct answer is 1/2. Here we seek an analogous ordering criterion for 
the bubbles. 
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neighborhood S of any point belonging to the "eternal fractal" image on Sq. Every such 
neighborhood will contain an infinite number of bubbles of all kinds and will be dominated 
by bubbles formed at very late times and having very small comoving sizes. The resulting 
values of Pj, obtained in the limit of bubble size e — ^ 0, will be the same in all neighborhoods, 
because of the universal asymptotic behavior of eternal inflation. 

It is clear that the same result will also hold in any finite-size region on Sq (provided that 
it contains at least one "eternal point"), and for any choice of the initial hypersurface Sq. 
Moreover, although we use the metric on Sq to compare the bubble sizes, the results are 
unaffected by smooth conformal transformations of the metric. Any such transformation 
will locally be seen as a linear transformation, which amounts to a constant rescaling of 
bubble sizes. In a sufficiently small neighborhood on Eg, all bubble sizes are rescaled in the 
same way, so the bubble counting should not be affected. 

The results obtained using this method are also independent of the initial conditions at 
the onset of eternal inflation.^ This is an attractive property, since the initial conditions are 
quickly forgotten in an eternally inflating universe.^ 

The calculation of bubble abundances, defined in this way, can be reduced to an eigenvalue 



problem for a matrix constructed out of the transition rates between different vacua 1^7 
This prescription has been tried in some simple models and appears to give reasonable results 



|24| . For example, if there is a single false vacuum, which can decay into a number of 



vacua with nucleation rates Tj , one finds 



Pj « r^-, (5) 

as intuitively expected. 



^ I assume that any vacuum is accessible through bubble nucleation from any other vacuum. Alternatively, 
if the landscape splits into several disconnected domains which cannot be accessed from one another, each 
domain will be characterized by an independent probability distribution, and our discussion will still be 
applicable to any of these domains. 

® An earlier suggestion of Ref. [19j was to define pj as the probability for a comoving observer to end up in 
a pocket of type j. jZbis probability, however, does depend on the initial state. 

^ The calculation in [ijj assumes that the divergence of geodesies is everywhere determined by the local 
vacuum energy density. This is somewhat inaccurate, since it ignores the brief transition periods following 
the bubble crossings and the focusing effect of the domain walls. The accuracyof the method is expected 
to be up to factors 0(1). A more detailed discussion will be given elsewhere |23| . 
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B. An equivalent proposal 

An alternative prescription for pj has been suggested by Easther, Lim and Martin [2^. 
They randomly select a large number of worldlines out of a congruence of comoving 
geodesies and define pj as being proportional to the number of bubbles of type j intersected 
by at least one of these worldlines in the limit of — ^> 00. 

As the number of worldlines is increased, the average comoving distance e between them 
gets smaller, so most bubbles of comoving size larger than e are counted. In the limit of 
00, we have e — * 0, and it is not difficult to see that this prescription is equivalent to 
the one described in the preceding subsection. (For a rigorous proof, see Note added ui jlS^.) 



C. Reference length Rj 

Our next task is to define the comoving reference scale Rj. The first thing that comes to 
mind is to set Rj to be the same for all bubbles. However, this is not enough. The expansion 
rate is different in different bubbles, so the physical length scales corresponding to Rj will 
not stay equal, even if they were equal at some moment. We could specify the times tj at 
which Rj are set to be equal, but any such choice would be subject to the criticism of being 
arbitrary. 

A possible way around this difficulty was proposed in At early times after nucleation, 
the dynamics of open FRW universes inside bubbles is dominated by curvature, with the 
scale factor given by 

aj(t)^t (6) 

for all types of bubbles. For example, for a quasi-de Sitter bubble interior, 

aj{t) ^ HTHin\i{Hjt), (7) 

where Hj = (SirCpj/Sy^^ is the expansion rate corresponding to the local vacuum energy 
density pj. The specific form of the scale factor at late times is not important for our 
argument. The point is that for t <C Hj'^ all bubbles are nearly identical, with the scale 
factor ©. (This is basically a consequence of the universal spacetime structure in the 
vicinity of the light cone t = 0.) 

n 

The proposal of [18|] is that the reference scales should be chosen so that Rj are the same 
at some small t = r. The choice of r is unimportant, as long as r <C Hj'^ for all j. Then, 



up to a constant, the physical length corresponding to Rj is 

Rj{t)=a^{t). (8) 

For times t 3> Hj^^ this can be expressed as 

R,{t)^H-'Z,{t), (9) 

where Zj is the expansion factor since the onset of the inflationary expansion inside the 
bubble [t ~ HJ^). 

Alternatively, Rj in (jHl) can be identified as the curvature scale. It is the characteristic 
large-scale curvature radius of the bubble universe. This definition makes no reference to 
early times close to the bubble nucleation: the curvature radius can be found at any time. 
It is, in principle, a measurable quantity. 

D. Continuous variables 

The prescription ^ can be straightforwardly generalized to the case when, in addition 
to bubbles, there are some continuously varying fields X: 

PP{X)^p,P,{X)R]{X). (10) 

Here, Pj{X) is the normalized distribution for X in a bubble of type j at t = r <C Hj^, 

P,{X)dX = 1. (11) 



This distribution can be calculated analytically or numerically, using the methods of 
Refs. 



V. DISCUSSION 



The above definition of the measure is just a proposal. We have not derived it from first 
principles. In fact, there is no guarantee that there is some unique measure that can be used 
for making predictions in the multiverse. How, then, can we ever know that we made the 
right choice out of all possible options? 

What I find encouraging is that even a single definition of measure that satisfies some basic 
requirements proved very difficult to find. In addition to being mathematically consistent, 
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we require that the measure should not depend on any arbitrary choices, such as the choice 
of gauge or of a spacehke hypersurface, and that it should be independent of the initial 
conditions at the onset of inflation. It would be interesting to know how much freedom is 
left by these requirements. In other words, how uniquely do they specify the proposal of 



Ref. |18|? 

The measure can further be tested by working out simple models and checking whether 
or not the resulting distributions are reasonable. So far, the proposal of Q| seems to have 
passed this test. 

A possible alternative, advocated in [l^, is to assume that "all infinities are equal" and 

set 

Pj = const (12) 

for all j. This prescription is clearly gauge-invariant and is independent of the initial con- 
ditions. The resulting measure, however, appears to be counter-intuitive. Intuitively, one 
expects that, everything else being equal, the probability assigned to a certain type of bubble 



3, 



should be proportional to the bubble nucleation rate. This is satisfied for the proposal of 
but if (fT^ is adopted, the probability would be completely independent of the nucleation 
rate. 

The ultimate test of any proposed measure will be a comparison of its predictions with 
observations. The first atteinpts in this direction have already produced some intriguing 
2^ I2S , 22 , 3, l^^, ^ QI , 13 • It seems safe to predict that we will hear more on this 



results 



subject in the future. 

I am grateful to Leonard Susskind for pointing out an error in the original version of 
this paper and for stimulating correspondence. I am also grateful to Jaume Garriga, Delia 
Schwartz-Perlov and Serge Winitzki for discussions and useful comments. This work was 
supported in part by the National Science Foundation. 
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